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^ ' Abstract 

£NJ ■ We show the instability of two self-tuning brane world models with gauge invariant 

linear perturbation theory. This general method confirms a known instability of the 
original model with vanishing bulk potential. We also show the dynamical instability 
' of a recently proposed self-tuning "smooth" brane model and its limit, the Randall- 

Sundrum model. Astonishingly, we also find instability under purely gravitational 

^ | perturbations. 

o ■ 

1 Introduction 

Superstring theories and especially M-theory have revived the idea that the universe is a 
^ ■ 3 + 1 dimensional hyper-surface, a 3-brane, of a 10 or 11 dimensional spacetime. This 
>• . fundamental spacetime could be a product of a four dimensional Lorentz manifold with an n 
dimensional compact space of volume V n . Then, the effective four dimensional Planck mass, 
M cS = y/l/(%irG N ) ~ 2.4 x 10 18 GeV is related by 

Ml = M? +2 V n . (1) 

to the 4 + n dimensional fundamental Planck mass, M { . Since M { is much smaller than 
M eff and may even be close to the electro-weak scale the long-standing hierarchy problem 
might be resolved. Clearly, this is an intriguing idea and it has attracted a lot of attention. 

Instead of a factorizable space-time with a compact extra-dimensional space, Randall and 
Sundrum |§](RS, see also (0)), proposed a five dimensional model, in which the metric on 
the 3-brane is multiplied by an exponentially decreasing 'warp' factor such that transverse 
lengths become small already at short distances along the fifth dimension. The brane is 
embedded in an Anti-de Sitter space, but unfortunately, a fine-tuning relation between the 
brane tension and the negative cosmological constant in the bulk has to be satisfied. 

Another possibility was discovered by Arkani-Hamed, Dimopuolos, Kaloper and Sundrum 
and by Kachru, Schulz and Silverstein || (ADKS-KSS). In this model, there is a 4D- 
Poincare invariant solution for each choice of the input parameters which is stable with 
respect to radiative corrections of the brane tension. This feature is called "self-tuning" and 
would explain why we observe a vanishing cosmological constant, i.e. it would solve the 
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cosmological constant problem ||. However, the model has singularities in the bulk which 
seem to be inevitable in self-tuning models of this kind ||10|| . It has also turned out that the 
model needs a fine tuning for the regularization of the singularity in the bulk [11]. A variant 
of these ideas was recently given by Kehagias and Tamvakis who introduce a self-tuning 
'smooth' brane model without singularities fl2"f , fl~3|[ (KT). 

These models are static. This immediately raises the question whether their essential 
features persist in time. For the ADKS-KSS model this is not the case: Arbitrarily close to 
the static self-tuning solution there are time-dependent solutions, which lead to an expansion 
or a contraction of the brane world; the model is therefore dynamically unstable M. 

Using gauge invariant linear perturbation theory we show in this paper that the ADKS- 
KSS model and the KT model are dynamically unstable against small variations in the scalar 
field and even homogeneous, purely gravitational perturbations. The instability of the KT 
model also implies the instability of the RS model because it is a certain 'brane limit' of 
the KT model. This limit does not modify the purely gravitational mode; therefore the RS 
model is also unstable under this mode. This result was also found in [ 16| . 

In section |2| we write down the dynamical equations for the gravitational and the scalar 
field which couples to gravity and to the brane tension. We derive linear perturbation 
equations which describe the evolution of small deviations from a static solution. In section 
[3] we give the perturbation equations in gauge invariant form and prove the instability of 
the KT model and the ADKS-KSS model. In the last section we present our results and 
conclusions. 



2 The Linearized Einstein Equations 

We consider a general five- dimensional action^ with a bulk scalar field <p with potential U(4>) 
and one brane with brane tension A (</>), 

S = [ d 5 x^\ (^R - Udcj)) 2 - U(<j>)) - [d 4 x v^A(0 b ) , (2) 



2k 2 

where g§ and #4 are, respectively, the determinants of the 5D metric jmjv and the 4D metric 
induced on the brane, g^, and denotes the value of the scalar field on the brane. The 
tension of the brane couples to the bulk scalar field in a way defined by the function A(0). 
The Einstein and scalar field equations are 

Gmn = k 2 (d M <pd N <P - g MN {^(d<pf + U(cj>)\ \ - ^& 9flv 6^ N 5(y) , 




1 dX y/\gi\ dU(<p) 

6 (y) = En • ( 3 ) 



2k~ u< > \g 5 



1 Conventions and Notation: We use a mostly positive Lorentzian signature ( h . . . +) and the 

definition of the curvature in terms of the metric is such that a Euclidean sphere has positive curvature. 
We write the 5D Einstein equations as Gmn = k 2 Tmn with the constant k related to the 5D Planck mass 
K 2 = I/4M5. Bulk indices are denoted by capital Latin indices and brane indices by Greek indices: are 
coordinates on the brane (t is the time coordinate and x 1 the spatial ones), and y is the coordinate along 
the fifth dimension such that the brane is located at y = 0. 
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To investigate the dynamical behaviour of brane world models, we use the following ansatz 
for the metric: 

9 = e 2A ^r } ^dx v + b(t, yfdy 2 . (4) 

Inserting this into the metric, one finds in the bulk (y ^ 0) the Einstein equations ( ' denotes 
d y and ' denotes d t ) 

3e- 2A (A 2 b 2 + Abb)-6A' 2 - 3A" + 3^- 

2 I 1 121,2-2.4 , L i/2 i 1,27 



= ^ + / + « (5) 

(-2i6 2 - A 2 b 2 - Abb - bb)e- 2A + 6A' 2 + 3A" - 3A' h - 

= ^(l^e-^-^-b'U^)) (6) 

-3&V 2A (i + i 2 ) + 6A /2 = k 2 Q0' 2 + ^0 2 - b 2 U(^ (7) 

-3i' + 3A^ = k 2 ^' , (8) 
and the scalar field equation 

e 2 ^ - 60 - 2A60) + - ^ + <f = • (9) 

V / 6 a0 

The 5(y)-terms in the action lead to jump conditions for the derivatives of the fields at the 
brane 

m,y)?a- = -^f 2 (10) 

m ,yy& = ^^f^- (ID 

To linearize these equations around a static solution (4>o(y), A (y), and 6 = 1)0, we set 

0(t,y) = My) + W(t,y) 

A(t,y) = A (y) + SA(t,y) and b{t, y) = 1 + Sb(t, y) . (12) 

The static solution has to satisfy the zeroth order field equations in the bulk, which reduce 
to 

6A' 2 = K 2 I(0') 2 -W), (13) 
-3A"-6A' 2 = Atf) 2 + U((f>) , (14) 
4A'0' + 0" = (15) 



2 We demand b(y,t) = constant, so with a redefinition of the y-coordinate we get 6=1. 
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and the zeroth order jump conditions at the brane 

KJS- = (16) 
m y =o- - 2k 2 d(j) • 

Many solutions of these equations are known; see for a general solving procedure for 
arbitrary U{<p). 

Now if we keep terms up to first order in 5, the (00), (ii), (44) and (04) components of 
the Einstein equations are 

-12A' 5A-3 5A" + 3A' 5b' 

= k 2 (<j)' 5<j)' + + 2U{<j) Q )8b 



-2 8 A - 5b)e- 2Ao + 12 A' 8A + 3SA" -3A' 8b' 

] -5<f)-2U((j) )5b) (19) 



2 I ,1 r-/ dU((j) 0j 



-3e- 2Ao 5A + 12A' M' 

k 2 Ufitf - ^i^H ~ 2U(<p )6bj (20) 
-35A+3A' 5b = K 2 (f)'5(f), (21) 



and the linearized scalar field equation is 

+^^64 > -4A' 8^ + ( p' 8b' . (22) 
At the brane the perturbation amplitudes have to fulfill the first order jump conditions 
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ny=0+ 1 d 2 X(<j)Q) 

ly =°- 2k 2 dd) 2 



(24) 



3 Perturbation Calculations 

The ideas and methods in this section come from cosmological perturbation theory (for 



an introduction, see ||14j|), originally developed to describe the formation and evolution of 
density fluctuations; here we use them to study questions of stability. Similar methods were 
used in [|lf| to derive the Friedmann equations of a brane world model. 

General relativity is invariant under diffeomorphisms (coordinate transformations). This 
means that for an arbitrary diffeomorphism (p and its pullback ip*, the metrics g and (f*{g) 
describe the same geometry. We look at small perturbations Sg around a static background 



4 



metric go and allow only diffeomorphisms which leave the background metric invariant. 
These 'infinitesimal' transformations can be described as the flow of an 'infinitesimal' vector 
field X and the perturbation of an arbitrary tensor field Q = Qo + SQ obeys the gauge 
transformation law 

SQ^SQ + L X Q , (25) 

where Lx denotes the Lie derivative in the direction of X. Taking the vector field to be of 
'scalar' type, i.e. 

X = X^ + X 4 «9 4 , X» = jf v i, v , (26) 

with an arbitrary scalar field £(x M ) and noting that the most general scalaif] perturbations 

M denotes ^ ) 



of the metric can be written as ( u denotes " 



Sg = 2e 2Ao {5 A + 8E >llv ) dx»dx v + 25bdy 2 + 2 SB >fl dx"dy, (27) 

we arrive at 

5 A -> SA + A' X 4 , (28) 

Sb -> Sb + X^, (29) 
6B jlt -> 5B^ + Xl + e 2A % M 

^SB -> 5B + X 4 + e 2Ao £ A , (30) 

,fJ,U 

^SE -> + (31) 

In brane world models there is a jump in A' Q at the brane, therefore (p8|) implies X 4 (x At , y = 
0) = and restricted to the brane, is gauge invariant. 

The scalar field perturbation 5(j) is defined by <p(t, y) = <po(y) + S(j)(t, y); according to ( p5|) 
one has 

50 -> 50 + <p' X 4 , (32) 
and again the restriction to the brane, S(pb, is gauge invariant. 

The gauge invariant terms for the perturbation amplitudes then are 

5A L := 5 A + A' (e 2Ao 6E, 4 - SB) (33) 
Sb L := 5b + (e 2Ao SE A - SB) 4 (34) 
50 L : = S(f) + (f)' (e 2Ao SE A -SB) . (35) 

Because of ( fflfy and fl3"T|), it is always possible to find a gauge with SB = SE = (longitudinal 
gauge). In this gauge, the metric is diagonal and has the form (f|), therefore it is much simpler 
to calculate the first order Einstein equations in the longitudinal gauge (|18|)-(|21"|) and later 
rewrite them in a gauge invariant form. In the longitudinal gauge each of the gauge invariant 
amplitudes above reduces to the corresponding gauge dependent term, for example SAl = SA 
etc. 



3 We only need perturbations of the scalar type because in T^ v we have just a scalar perturbation of the 
scalar field <f>. 

4 The first two terms, SAl and Sb^, correspond to the Bardeen potentials in cosmological perturbation 
theory if we identify our extra dimension with the time dimension and our 4D Minkowski space with the 3D 
space of constant curvature. 
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Two other gauge invariant combinations of the perturbation amplitudes will be useful in 



the following: 

: = SA' + A' 5b - — (f>' 6<p (36) 
$ := ^d^-M' + ^db, (37) 
where the gauge invariance of $ is a consequence of the zeroth order equation 3Aq = — k 2 0' o 2 , 



which follows from (|T3"D and (fH]). In the definitions ( |36"D and (|37| ) we could of course equiv- 
alently use the gauge invariant amplitudes SA L , 56^ and 50^; the SB and terms cancel. 

Using the zeroth order equations (p~3|)- (p~5|) we can absorb the terms with the potential 
U((f>o) in the first order equations into combinations of the gauge invariant terms $ and $. 
In gauge invariant form, the (00), (00) + (ii), (44) and (04) components of the first order 
Einstein equations (@)~(E3) are 

(38) 
(39) 

(40) 

, (41) 

and (|2^), written in gauge invariant form, gives 

64," L = Sh + #;* + (2^^ - 4^) Sb L 

+ (^^ + ^)«L-4^% + ««i. (42) 

The first order jump conditions fl23l) and (pl|) are already gauge invariant because of X 4 (x M , y = 
0) = and X 4 G C 1 . 

The procedure to find suitable perturbation amplitudes 5Al, 5bi and 5(pL is to solve 
(|38|)-(|4"2|) in the bulk, i.e. for y < and y > 0, and then glue the solutions continuously 
together at y — and check whether they satisfy the jump conditions. 



AA'^ + V = 

-3e~ 2Ao 5A L - 12A ^ + k 2 $ = 

^ = 



3.1 A Solution with 5bL = 

Now we derive a solution to the bulk equations which we will later use to prove the instability 
of the KT and the ADKS-KSS self-tuning model. Equations ( |38"D and fl4l|) imply 



V(t,y) = K ie - 4Ao{y) , (43) 

where K\ is an integration constant. We try to find perturbations linear in t, and to simplify 
matters we set SbL = 0. (f4"0"|) then yields 

$ = ^-A' V = - 2 A K ie - iA ^ . (44) 



According to the definitions of $ and \I> (P6"|)-(|TF|), this implies 

12 A t 

~2 
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The subsequent integration gives 

/A A' 
-£e- 4A °dy + k 2 (t)<f>' , (46) 

with another (time dependent) integration constant k 2 . From the definition of \& we can find 
5A' L and then the perturbation amplitude SAl must be 

SA L = KtJ 0' o 2 ^ ^e" 4Ao dy^j dy + A' k 2 (t) ~ Ki J ^ M dy + k 3 {t) . (47) 

One can check that the scalar field equation is also satisfied. 

Since the solutions were assumed to be linear in t, the integration constants k 2 (t) and 
ks(t) have the form ki(t) = + hit, and we get a five-parameter family of perturbation 
amplitudes 

SA L = K,J 0' o 2 (| ^e-^dy^dy 

+A' K 2 + A' h 2 t -KiJ e~ 4Ao dy + K 3 + h 3 t (48) 



/A A' 
~^e~ iAo dy + <P',K 2 + <t>> Q h 2 t. 
ro 



(49) 



We are only interested in time-dependent perturbations, so we set Ki = and obtain 

SA L = A' h 2 t + h 3 t (50) 
5<Pl = <P' h 2 t. (51) 

Note that for hi = all the gauge invariant terms are zero, and the solution is equivalent 
to the unperturbed, static solution A = A , b = 1, = O an d SB = SE = 0. 

3.2 Instabilities in the KT-Self-Tuning-Model 

In the self-tuning model proposed by Kehagias and Tamvakis there is no brane term propor- 
tional to S(y) in the five dimensional action (^), but a 'smooth' brane is dynamically formed 
by the scalar field coupled to the Standard Model. So we have to set A(0) =0 to study the 
stability of this model. Now the field equations are valid everywhere and there are no jumps 
in the derivatives A(t, y) and <j)(t, y), the jump conditions are trivial. The potential is chosen 
to be 

U(<P) = V (1 - k 2 sin 2 (/30)) with k 2 = 1 + ^ . (52) 

The static solution <po(y) has the kink-like form 

2 

P 

where the parameters are related and restricted by 



o = — arctan ^tanh(cM//2) j , (53) 



2V 

Q2 " 2 ^ 2= 3M3(P-1) and fc2>1 ' (54) 
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and the static warp function A (y) is 



(1 — k 2 ) 

Mv) = 4 In cosh(ay) . (55) 
For this solution, the four-dimensional Planck mass has a finite value 

p ~ 5 7-oc y ~ <k 2 - 1) r ((P + l)/4) ' (56) 

and the geometry is nowhere singular for k 2 > 1. The solution depends on two parameters, 
k 2 and Vo. Standard model physics is represented by V , and for every value of V there is a 
solution with a flat 4D hyper-surface, i.e. the model displays self-tuning.f] 

The solutions (|5CiD , (|5l| ) of the perturbed bulk field equations directly show the dynamical 
instability of the KT-model since the bulk equations are valid everywhere and we have no 
jump conditions for this model. Now we will discuss two special perturbation amplitudes 
that form part of these solutions. 

First we consider a perturbation with = in (|50|) 



5Al = A' h2t and 5<pL = 0o^2^ , (58) 

i.e. the model is unstable against small perturbations in the scalar field. The instability 
could have the unpleasant effect that the 4D Planck mass becomes time dependent, but a 
second order analysis is needed to settle this point. To first order in 5 A we have 

/oo /*oo 
e 2A(t, y ) dy M 2 q +M I 2 5A(t, y)e 2My) dy , (59) 
-oo J — oo 



where M p o denotes the unperturbed 4D Planck mass from (p6|). The special perturbation 
yields 

/oo 
2A' (y)e 2Ao ^h 2 tdy = M 2 , (60) 
-oo 

and there is no change compared to the static solution in this approximation. It is possible 
that the (unknown) dynamical solution of the full field equations, which has (|58D as first 
order time expansion coefficient, also leads to usual 4D gravity on the brane. However it is 
not clear to us whether energy is still conserved on the brane in such a dynamic geometry. 

Now we set h 2 = in (|50"D. Then the perturbation is 

5A L = h 3 t and 5<p L = . (61) 



There is an interesting brane limit to the smooth solution (53) and (p5|), defined as Vo,a — > oo and 
k 2 — > 1 so that £ 2 := 8n 2 Vo(k 2 — l)/3 remains finite. In this limit, A$(y) — > — £|y|/4, the scalar field and 
potential part of the 5D action goes to 

1^2 , TT ,^\ ... f°° , 3 £ 



2 . ■m)dy^j_ oo ^ + ^s(y))^y- (57) 

Now it corresponds to the action of the Randall-Sundrum model [0, with bulk cosmological constant 
A = — 3k _2 ^ 2 /8 and brane tension Vb rane = 3k~ 2 S;/2. Note that the RS fine-tuning relations are fulfilled 
automatically! 
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Thus there is also an instability if the scalar field <p is n ot perturbed. The thick brane model 
is unstable even under homogeneous gravitational perturbations^. The perturbed 4D Planck 
mass is now time-dependent 

/oo 
2e 2Ao{y) ht dy = M p 2 (l + 2h 3 t) . (62) 
■oo 

The gravitational constant G is proportional to + 2h 3 t = 1 — h 3 t. This need not be 

in conflict with observation if /i 3 is small enough, but it surely is an unpleasant effect of the 
instability. Although the physical interpretation of these instabilities is not complete, we see 
that the KT self-tuning model requires a fine-tuning to be stable in time. 



3.3 The Instability of the ADKS-KSS Model 

In the ADKS-KSS self-tuning model (see JF| and ||) the bulk potential U((f>) is zero and the 
static solution is 

A (y) = llnfl-M^ (63) 



4 V y 



Ami -J* 1 



= ^+\/772 ln ( 1 --) (64) 



A(0) = ce" 2 ^/^ (65) 

y c := h- 2K ^ . (66) 

In the following, we will obtain the same instability with the more general method of gauge 
invariant linear perturbation theory. 

We first consider a perturbation with h 3 = 

SAl = A' Q h±t and 5<pL = (fi'oh±t . 

The field equations are only valid in the bulk, i.e. for y ^ 0. Therefore we actually have two 
integration constants, h + where y > 0, and h- where y > 0. A and 0o are Z 2 symmetric, 
A' and 0o are Z 2 antisymmetric with a jump at y = 0, consequently we have to set h : = 
h + = —h^ to get continuous perturbation amplitudes 5Al and 5<pL'- 

— ht 

5A L (t, y) = sgn(y)A' (y)ht = j^e~ 4A ^ (67) 

H L = m (y)<&(y)ht = —.fJie-iMy) . (68) 

y c V 4:K 2 



These perturbations also satisfy the first order jump conditions ( P3[ ) and fl2"4j), so they show 
the dynamical instability of this model against small changes in the scalar field. The solution 
is analogous to the dynamical solution found in [], the perturbations 5A L and S(f>L are 

6 In the brane limit this mode survives unchanged and shows that the Randall-Sundrum-model is unstable 
under purely gravitational modes. This result was also found in Q. 

7 The dynamical solutions in |j] are obtained with a diffeomorphism and orbifold projection from the 
static solution. Therefore it is possible to transform away the perturbation amplitudes locally in the bulk. 
The transformation, under which the metric remains diagonal, corresponds to the vector field with X A = 
—htsgn(y) and £ = ht J exp(—2Ao(y))sgn(y)dy. Note that this is not a gauge transformation because 
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proportional to the first order terms in a time expansion of their dynamical solution around 
(HI) and d64|) . See for the physical interpretation of this instability. 

The homogeneous gravitational mode (|JT]) is continuous and fulfills the first oder jump 
conditions since the jumps are zero and S^l is zero as well. Therefore this model is also 
unstable under homogeneous gravitational perturbations. 



4 Results and Conclusions 

Using linear perturbation theory we have demonstrated the instability of some brane models. 
In particular, we have reproduced the instability of the ADKS-KSS model J7J, |8| obtained 
by a diffeomorphism and orbifolding in ||. Linear perturbation theory also works for the 



smooth KT brane model [|13| and presumably also for brane world models with non-trivial 
bulk potential. 

The KT self-tuning model is not stable against small perturbations of the scalar field and 
it is even unstable under a purely gravitational mode since there is a metric perturbation 
linear in time which solves the field equations with the unperturbed scalar field. This mode 
causes the 4D effective Planck mass to be time dependent. The physical interpretation of 
these instabilities is not clear to us; it is possible that they affect energy conservation on the 
brane, as in the ADKS-KSS model. The instabilities show that the KT self-tuning model 
requires a fine-tuning to obtain the static solution. 
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